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This  paper  presents  a  framework  for  synchronously  coupling  wave,  curreni,  sediment 
transport,  and  seabed  morphology  for  the  accurate  simulation  of  multi-physics  coastal 
ocean  processes.  The  governing  equations,  which  represent  models  that  are  commonly 
adopted  in  practical  simulations,  are  discretized  using  finite-difference  methods.  The 
resulting  system  is  validated  against  analytical  solutions.  In  order  to  test  the  performance 
of  the  proposed  framework  and  the  numerical  methods,  dam-break  flow  over  a  mobile-bed 
and  evolution  of  a  wave-driven  sand  dune  are  simulated.  The  interactions  among  waves, 
currents,  and  seabed  morphology  are  illustrated. 

©  2008  Elsevier  B  V.  All  rights  reserved. 


1.  Introduction 

Nearshore  hydrodynamics  has  a  direct  impact  on  societal  issues  such  as  coastal  engineering,  environmental  protection, 
recreation,  and  military  operations  (e.g.,  [1,2]).  Wave,  current,  sediment  transport,  and  morphology  are  important  processes 
within  coastal  and  estuarine  settings,  and  consequently  a  number  of  models  have  been  developed  to  simulate  and  predict 
their  behaviors  in  the  past  few  decades.  For  example,  models  have  been  designed  to  forecast  global  currents,  salinity,  sea 
level,  temperature,  and  turbulence  distributions  (e.g.,  [3.4])  and  nearshore  wave-driven  currents  (e.g.,  [5]).  In  addition,  a 
number  of  models  have  been  designed  to  simulate  surface  wave  propagation  (e.g.,  [6,7]).  Models  have  also  been  proposed 
to  predict  sediment  transport  and  seabed  morphology  (e.g.,  [8,9]). 

Since  actual  coastal  ocean  flows  are  multi  physical  processes  with  complex  interactions  among  hydrodynamic,  sediment 
transport,  and  morphological  phenomena,  the  consideration  of  feedback  between  these  individual  physical  phenomena  is 
necessary  to  accurately  simulate  the  flows.  Model  coupling  has  been  implemented  for  the  prediction  of  many  actual  prob¬ 
lems;  for  instance,  a  wave  model  has  been  incorporated  into  a  primitive-equation  model  to  study  wave-current  interaction 
[  1 0).  This  study  demonstrated  that  decoupling  currents  and  waves  leads  to  inaccurate  water  elevation  prediction  Wave  and 
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morphology  models  have  been  added  into  the  Delft  Hydraulics  model  (DELFT3D)  (e.g.t  [9]).  Efforts  have  also  been  made  to 
couple  the  Simulating  WAves  Nearshore  (SWAN)  model  with  a  current  model  to  simulate  sediment  transport  (e.g.,  [11]). 
Nevertheless,  strictly  speaking,  until  now  the  coupling  between  currents,  waves,  sediment  transport,  and  morphology  has 
been  non-synchronic  or  merely  partially  implemented.  For  example,  the  wave  computation  is  not  simultaneously  carried 
out  with  that  for  currents  and  morphology  in  simulations  described  in  (9);  also,  in  coupling  of  SWAN  with  the  Advanced  Cir¬ 
culation  (ADC1RC)  model  [12],  the  wave  component  is  updated  only  once  in  solving  the  current  and  sediment  transport 
equations. 

In  order  to  understand  complex  phenomena  within  coastal  ocean  flows,  it  is  necessary  to  fully  couple  waves,  current,  and 
sediment  transport/seabed  morphology.  To  this  end,  a  number  of  investigations  have  been  made  on  numerical  methods  for 
such  coupling.  Rogers  et  al.  [  1 3]  designed  a  Godunov-type  method,  together  with  an  adaptive  grid  scheme,  to  simultaneously 
simulate  ray-type  wave-current  interaction.  Murillo  et  al.  [14]  combined  inviscid  shallow  water  and  solute  equations  into  a 
single  system  and  proposed  a  finite  volume  method  based  on  the  Roe  average.  Hudson  and  Sweby  [15,16]  made  a  system¬ 
atical  study  of  different  formulations  to  couple  currents  and  morphology.  Rosatti  and  Fraccarollo  [17]  constructed  a  Godunov 
scheme  for  currents  over  a  mobile-bed.  An  important  issue  is  discretization  of  the  source  terms  in  the  coupling  systems  (e.g., 
[18]).  A  so-called  C-property  is  proposed  to  enforce  balance  between  numerical  fluxes  and  non-homogeneous  source  terms 
in  cases  of  steady  flows  (e.g.,  [19]).  Nevertheless,  in  previous  investigations,  the  coupling  has  been  usually  between  two  or, 
occasionally,  three  components  among  wave,  current,  sediment  transport,  and  morphological  processes. 

In  order  to  accurately  take  their  interactions  into  account,  this  paper  presents  a  framework  to  synchronously  couple  wave, 
current,  sediment  transport,  and  morphological  processes.  This  coupled  system  consists  of  governing  equations  that  repro¬ 
duce  typical  models  used  for  both  research  and  engineering  problems.  The  wave  action  equation  is  used  for  wave  propaga¬ 
tion.  The  shallow  water  equations  are  employed  to  describe  currents  and  storm  surge.  Sediment  transport  is  described  by  a 
convection-diffusion  equation.  Morphology  evolution  is  depicted  by  the  Exner  equation.  In  this  paper,  all  of  these  equations 
are  discretized  using  common  schemes,  which  advance  in  time  synchronously,  or  with  same  time  steps.  The  resulting  system 
is  validated  using  analytical  solutions.  In  order  to  demonstrate  the  performance  of  the  framework,  simulations  are  presented 
of  dam-break  flow  over  a  mobile-bed  and  wave-driven  flow  over  a  sand  dune. 

The  paper  is  organized  as  follows.  Section  2  presents  the  coupled  system  with  governing  equations  for  waves,  currents, 
sediment  transport,  and  morphology,  together  with  a  discussion  on  their  coupling.  Section  3  deals  with  discretization  of  the 
system.  Section  4  discusses  validation  of  the  system  and  two  example  simulations.  Section  5  concludes  the  paper. 

2.  Governing  equations 


The  governing  equation  for  wave  action  is  given  in  a  conservation  form  as  (e.g.,  [20]) 
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where  f  is  the  time,  x  and  y  are  the  Cartesian  coordinates,  a  is  the  frequency,  and  0  is  the  angle  of  the  wave  propagation  direc¬ 
tion.  N  is  the  wave  action,  Q*  and  Qy  are  respectively  the  wave  speed  inx-  andy-direction  in  the  physical  space  (x,y),  C„  and 
C0  are  respectively  wave  speed  in  a-  and  0-direction  in  the  spectrum  space  ( a ,  0).  S  is  a  source  term  that  represents  the  com¬ 
bined  effects  of  wind  and  other  processes.  Using  linear  wave  theory,  it  has  been  shown  for  short  waves  that  (e.g.,  [20]) 
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and  for  long  waves  that 
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In  the  above,  U  and  V  are  the  current  velocities  in  x-  and  y-directions,  respectively,  H  is  the  water  depth,  and  g  is  the  grav¬ 
ity  constant.  For  convenience  in  computation,  only  short  and  long  wave  situations  are  considered  here.  Short  waves  are  de¬ 
fined  as  those  with  kH  >  n,  and  long  waves  as  those  with  kH  <  zr/1 0.  Here  k  is  the  wave  number.  For  further  details,  readers 
are  referred  to  [20,21]. 
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The  shallow  water  equations  consist  of  the  equation  of  mass  conservation  (e.g.,  [ 20] ) 
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where  vf  is  the  turbulence  eddy  viscosity,  and  p  is  the  density.  //  is  the  water  surface  elevation,  r*  and  xsy  are  the  surface  fric¬ 
tion  stresses  in  x-  and  y-directions,  respectively,  and  xbx  and  rb  are  the  bottom  friction  stresses  in  x-  and  y-directions,  respec¬ 
tively.  S**,  Sxy,  and  Syy  are  the  radiation  stresses  resulting  from  waves.  It  is  readily  seen  that,  for  the  case  of  a  flat  bottom,  Hdtjl 
0x«0(H2/2)/0x  and  Hdr]ldy  -  0(H2/2)/0y,  and  the  shallow  water  equations  will  be  in  conservation  forms  with  the  friction 
stresses  as  source  terms. 

The  turbulence  eddy  viscosity  is  calculated  as  follows  (e.g.,  [22,23]): 


v,  a,gn2H2  ^\JU2  4  V2. 


(5) 


where  ar  is  a  constant  ranging  from  0.3  to  1,  n  is  the  Manning  friction  coefficient.  It  is  proposed  to  treat  wave  effects  on  bot¬ 
tom  friction  as  follows  (e.g.,  [24]): 


tj  f>gn2H  ,/3U\U2  4  V2  +  \pfwU2w.  pgn2H~'/3V\j  U 2  4  V2\pfwul 


(6) 


where  Uw  -  An/Sin  h{kH),  the  near-bed  wave  orbit  velocity,  and  fw  is  the  wave  friction  factor.  Here,  A  is  the  wave  amplitude. 
In  the  case  of  short  waves,  kH  >  1  and  a  =  v  gk,  thus  Uw  «  0,  and 


rj  pgn2H  '/3U\U2  4  V2 


pgn2H 


,V</u2~i~V2. 


(7a) 


In  the  case  of  long  waves,  kH  <5C  1  and  a  =  kx  gH,  thus  Uw  ^  AxgjH.  Since  £  -  pgA2\ 2  (£  is  the  wave  density  and  it  equals 
N\a),  Uw  %  y/2 WfpoH.  As  a  result,  one  has 


pgn 


H  '/3U\  U2  4  V2  +fwN 


a  H  ’  Ty 


pgn2H  '  3Vy  U2  4  V2  4 


For  the  surface  stress,  a  widely  used  formula  is  (e.g.,  [25,26]) 


CoPairWxW 


CoPmrWyW 


(7b) 

(8a) 


where  Wx  and  Wv  are  wind  speed  in  x-  and  y-directions,  respectively,  W  is  the  total  wind  speed,  W  =  y  W2  f  W2,  paiT  is  the 
air  density,  and  the  drag  coefficient  is  given  as 


Co 


1.2875  x  10  \  W  ^  7.5  m/s, 

(0  8  4  0.065W)  x  10  3.  W  >  7  5  m/s. 

Assuming  linear  wave  theory,  the  radiation  stresses  for  short  waves  are  given  as  (e.g.,  [20]) 
g 


Jxy  ^yx 

g 

-yy  2 
and  for  long  waves, 


~  Jj  N  cos2  0  dado 

|  JJ  NsmOcosO  dadO , 


Jj  N sin2  0  dadO , 


g  JJ  N( cos2  4  0.5)  dadtf, 
[  II  N  sir 


f  sin  0  cos  0 


Syy  g  j j  N{ sin2  0+05)  drrdd. 


(8b) 

(9a) 

(9b) 

(9c) 

(10a) 

(10b) 

(10c) 


The  sediment  transport  equation  is  a  convection-diffusion  equation  with  a  source  term,  for  which  the  conservation  form 
is  (e.g.,  [22]) 


0HC  dHUC  dHVC 


dt 


+ 


dx 


+ 


Dc 


d2HC  4  d2HC 
dx2  dy 2 


-  0L(O{C  -  C) 


(11) 
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where  C  is  the  concentration  for  sediment,  C*  is  the  sediment  concentration  under  equilibrium  conditions,  Dt  is  the  diffusivity 
coefficient,  a  is  the  coefficient  of  suspended  load,  and  (o  is  the  settling  velocity.  The  seabed  morphology  is  controlled  by  the 
Hxner  equation  (e.g.,  [27,28]): 

SHi,  +  V2)*  ,  cW(U2  +  V2)9  Q 

cT  ^  dx 

Here  Hb  is  the  elevation  of  the  seabed,  and  R  and  q  are  constants.  Typically,  B  depends  on  flow  velocity,  water  depth,  sediment 
grain  sizes,  and  other  factors,  and  q  is  usually  in  the  range  of  0.5  ^  q  ^  1 .5,  t]  *  H  +  Hb.  In  this  paper,  q  -  1  is  used. 

Eqs.  (1),  (4),  (11)  and  (12)  comprise  a  coupled  non-homogeneous  system  of  conservation  laws,  which  reproduces  the 
framework  of  well-known  models.  Eqs.  (1)  and  (4)  can  be  viewed  as  a  result  of  flow  decomposition  in  the  Navier-Stokes 
equations  (e.g.,  [20]).  The  wave  action  Eq.  (1)  is  employed  in  models  such  as  in  SWAN  [6].  The  shallow  water  Eqs.  (4)  are 
used  for  models  such  as  SHORECIRC  [5].  The  sediment  transport  Eq.  ( 1 1 )  and  morphology  evolution  Eq.  ( 1 2)  are  also  widely 
employed  in  engineering  (e.g.,  [22,29]). 

The  interactions  among  wave,  current,  sediment  transport,  and  morphology  can  be  seen  in  the  system  equations;  in  Eqs. 
(2)  the  wave  field  is  coupled  with  the  current  through  the  velocity,  Uand  V.  As  indicated  in  Eqs.  (4)  the  current  is  affected  by 
the  wave  field  through  the  radiation  and  bottom  stresses,  Sxy,  Syy  and  tj,  xby,  respectively,  and  by  sediment  transport  and 
morphology  evolution  through  turbulence  eddy  viscosity  vr  and  bottom  elevation  Hb.  Both  sediment  transport  and  morphol¬ 
ogy  are  directly  related  to  the  current  through  the  velocity  field  as  shown  in  Eqs.  (11)  and  (12).  It  should  be  noted  that,  as 
seen  in  Eqs.  (2)  because  wave  action  occurs  in  spectral  space  as  well  as  physical  space,  each  term  on  the  right  hand  side  of  Eq. 
( 1 )  is  coupled  with  the  current  through  the  velocity  field.  Furthermore,  sediment  concentration  and  morphology  are  directly 
related  to  each  other  (e.g.,  [22]).  In  this  paper,  however,  in  view  that  their  governing  Eqs.  (11)  and  (12)  have  different  un¬ 
knowns,  which  are  sediment  concentration  and  bed  elevation,  sediment  transport  and  morphology  are  indirectly  coupled 
through  the  velocity  field. 

3.  Numerical  methods 


The  coupled  system  of  waves,  currents,  sediment  transport,  and  morphology  are  discretized  as  follows.  For  wave  action 
Eq.  (1),  the  following  first-order  accurate  extension  of  the  one-dimensional  Lax-Friedrich  scheme  [30]  is  used: 


r/.fc'.r  i.i 

i  At( 


<(*)'  . 


(13) 


where  At  is  the  time  step,  Dm  is  the  central  difference  operator  with  respect  to  the  m-direction  (m  -  For  instance, 

DDijjcj  -  (()i*ij>r()i  i  j.k.i)/(2Ax),  Ax  being  the  grid  spacing  in  the  x-direction.  The  shallow  water  Eqs.  (4)  and  the  sediment 
transport  Eq.  (11)  are  solved  using  the  following  MacCormack  scheme  [31  ]: 


Hy  Hi;  -  AfD,  (HU^J  -  At D,  (HV)”y 

(HU)lj  (HU)"j  -  AtD,  (HU%  AtDj  ( HUV )■  -  gAtH,"  1/2jD,  (H  i  Hb)?j  +  2AfD„  ((v,H)f  v2jD,  U”j) 


+  AfDJ+((v, H)^1/2Dj  l/y)  4  AtDj((v,H)"jD,y"j)  + 
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(W)y  (HV)1j  -  AtD ,  ( HUV )”j  -  AtD).(HV2)1j  -  gAtH"j_,  2D,  ( H  +  Hb)"j  f  AtDj((v,H)yD2l/y) 


AtD,  ((V,H)y  ,  2D,  V"j)  +  2AtDj+((v,H)y  ,  2d;  V£)  4 


At 


At, 


d:st  -  —  d;s;  , 

p  1  p  j  w ij 
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H?;’  j(Hy  *  «y)  -  yDl+(HU)y 

(ml"  \((HU)l  f  mi)-  y  Dn(HU%-^DJ+(HUV)l-^Hl,2jDl+(H  t  -  H*)y 


A tDh((v,H);  l/2JD,  U,v)  +yDJ+((V,H)y..1/2D;  (Ty)  4  y ^((v.H^D.yv) 
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Af(tJ'  rj'.)  At 
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(14a) 


(14b) 


(14c) 


(14d) 
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(HV) 


n  +  1 

ij 


j((HV)y  +  (HV)l)  -  ^-D,+  (HUV)‘j-YDb(HV%  ~^H]jUj2Dh(H  4  H»)y 
+  fDi((v,H);jDjU;j)  +yDlt((v,H);  i/2jD,  Vy)  +  AtD,+  ((V,H)y_1/2Dj  Vy) 


At(T^  -  gt j 
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D,S- 


yy<p 


(14f) 


(HC)V  (HC)"j  -  AtD,  (HUC)"j  -  AtD,  (HVC%  +  AtD14  (d,  (HC)£)  +  AfDJt  (D,  (HC)£),  (15a) 

(HC)f;1  1((HC)J1  +  (HQIj) -^■Du(HUC)\i-^Dh(HVC)\i ryD14(ft  (HC)V)  -  y  DJ4(D,  (HC)V).  (15b) 

In  the  above,  Dm  (),j  and  Dm+(),j  (m  -  ij,k,  /)  are  backward  and  forward  differences,  respectively.  The  values  at  half  nodes  are 
evaluated  using  an  arithmetic  average.  For  example,  0j+i/2j  ■  ((),j  +  ()J>iJ)/2.  The  morphology  Eq.  (12)  is  discretized  using  the 
following  scheme: 

H%'  =  BAtD,[L/(U2  4  V2)]'’  -  BAtD,[V(U2  4  V2)]^.  (16) 

The  time  steps  for  integrating  the  wave,  current,  and  sediment  transport/morphology  equations  are,  respectively 


N  (kg/s1) 


o(1/S) 


Fig.  1.  Numerical  solution  of  wave  propagation  given  by  (17a)  and  (17b).  T  -  3.5  s.  Error  j,N  N°  ^-  3.37  x  10  2.  where  N  and  N°  are  respectively 
numerical  and  analytical  solution. 
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CFL  •  Ax  CFL  Ay  CFL  ■  A  a  CFL  •  Afl] 

^-gXijkl  Cgyijjlj 

CFL  Ax  CFL  Ay  von  min(Ax2  Ay2) 


mm 

U 


\Utj\  +  V^  l^l  +  y/gftj' 
CFL  Ax  CFL  Ay) 


m  ’  IVijI 

where  CFL  and  von  are  constants.  The  time  step  for  the  coupled  system  is 
At  mm{(At)1t(At)2,(Af)3} 


(17a) 

(17b) 

(17c) 


( 1 7  d ) 


4.  Numerical  examples 

4.1.  Validation 

The  numerical  methods  presented  in  the  previous  section  and  the  corresponding  computer  code  are  validated  using 
various  test  problems  with  analytical  solutions.  The  first  test  problem  is  long  wave  propagation  in  a  uniform  current  with 


Fig.  2.  Solution  for  surge  problem  Eqs.  (18a)  and  (18b).  T  -  2  s.  Lines  -  exact  solutions,  circles  -  numerical  solutions. 
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U  1  m/s,  V  1  m/s,  H  1  m  (18a) 

and  the  source  term  given  as 

S  4asN((x  -  t)2  +  (y  t)2  t2)(x  t  y  -  2t  a,t  C^x  -  t)  -  Cgy(y  -  t)) 


t  a6N 


(t. 


1  1  n 

- - - -  +  —  2a2 (a  -  a3) 

exp  -a6t  <? 


Q4NC0  cos  a4f 
2  +  sina4d 


(18b) 


A  solution  of  Eq.  ( 1 )  is 

N  a7(  1  exp{  a6t})(2  f  sin a40)  exp{  a2(a  -  a3)2}  exp{-a5((x  -  t)2  f  (y  -  t)2  -  a^2)2},  (18c) 

which  is  a  moving  circular  ring  with  a  diameter  increasing  with  time.  Let  a,  *  1,  a2  *  001,  a3  *  2.5,  a4  -  1,  a5  «  0.01,  a6  -  0.05, 
a7  -  1,  Ax,  Ay  -  0.4  m,  and  CFL  -  0.5.  The  numerical  solution  is  plotted  in  Fig.  1  with  the  error  given  in  the  figure  caption. 

The  second  test  problem  is  a  one-dimensional  surge-current  without  viscosity,  surface,  and  bottom  friction.  The  initial 
condition  is 


x<0.  U  2  m/s  V  0  m/s,  H  5  m  (19a) 

x>0  U  2  m/s,  V  0  m/s,  H  3  m  (19b) 

Letting  Ax  -  Ay  -  0.2m  and  CFL  -  0.95,  the  numerical  solution  is  obtained  as  shown  in  Fig.  2,  together  with  the  corresponding 
exact  solution  (e.g.,  [32]). 


X  (m) 


Fig.  3.  (a)  Numerical  solution  for  sediment  transport  from  a  point  source.  T  -  4  s.  Error  C  C°  x  -  1.35  x  10  2.  where  C and  C°  are  respectively  numerical 
and  analytical  solution,  (b)  Morphology  evolution.  T  -  5  s.  Lines  -  exact  solutions,  circles  -  numerical  solutions. 


2942 


H.S.  Tang  et  al./Cammun  Nonlinear  Sci  Numer  Simulat  14  (2009)  2935-2947 


The  third  test  problem  is  an  initial  value  problem  for  sediment  transport  with  a  point  source  that  has  the  following  solu- 
tion  ((33 ]): 


K  (x  —  Ut)2  (x  -  Vt)2 

4nDct CXP  4Dct  ~  4Dct 


(20) 


The  point  source  is  initially  located  at  (x,y)-(0,0),  and  K-  50,  U,  V -  1  m/s,  and  Dc-  1  m2/s.  With  Ax,  Ay -0.4  m, 
CFL  -  0.95,  and  von  -  0.2,  the  numerical  solution  is  obtained  as  shown  in  Fig.  3a.  The  last  test  problem  is  an  initial  value  prob¬ 
lem  of  one-dimensional  seabed  evolution  from  initial  elevation 
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Fig.  4.  Solutions  for  water  surface  and  bed  morphology  elevations  (z)  of  dam-break  flow  aver  mobile-bed. 
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H„  =  l +cos(^),  (21a) 

with  BU{U2  +  V2)q  -  1/H  in  Eq.  (12)  (e.g.,  [ 29)).  An  exact  solution  is  (e.g.,  [27,29]) 

H‘  i+co^tc(','<^«1))'  (2lbl 

The  numerical  solution  is  presented  in  Fig.  3b  with  Ax  -  0.2  m,  and  CFL-0.95.  All  of  these  tests  indicate  that  the  numerical 
methods  and  the  computer  code  work  correctly. 
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Fig.  5.  Solution  of  wave  action  for  the  wave-driven  dune  problem. 
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4.2.  Applications 

The  first  application  represents  the  interaction  between  surge  and  bed  morphology  in  an  idealized  dam  break  experiment 
over  a  movable  and  initially  flat  bed,  as  described  by  Spinewine  and  Zech  [  34],  A  horizontal  glass-walled  flume  of  rectangular 
cross  sectional  geometry  was  used.  The  initial  bed  profile  featured  flat  horizontal  reaches  upstream  and  downstream  of  a 
gate,  with  a  saturated  bed  of  constant  thickness.  In  the  upstream  reach,  a  water  layer  of  0.1  m  was  retained  and  released 

at  time  0  s  by  raising  the  gate  rapidly.  The  problem  is  formulated  as  the  following  initial  value  problem: 

x  ^  0  U  =  0  m/s,  V  =  0  m/s,  H  =  0.1  m,  Hb  =  0  m  (22a) 

x  >  0  :  U  0  m/s.  V  0  m/s,  H  0.002  m  Hb  0  m.  (22b) 


X  (m) 


Fig.  6.  Solution  of  current  for  the  wave-driven  dune  problem. 
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Here  a  thin  layer  of  water  is  given  at  x  >  0  for  convenience  of  computation  and  the  Lax-Friedrich  scheme  [  30]  is  used  for  the 
surge  with  the  following  parameters:  n  -  0.025,  B  «  0. 1  1  /s,  Ax  *  0.05  m,  CFL  -  0.95,  and  von  -  0.5. 

The  numerical  solution  is  obtained  as  shown  in  Fig.  4.  In  the  figure,  H0  -  0.1  m,  t0  =  v  H0/g,  and  numerical  solutions  and 
experiment  measurements  for  both  water  surface  and  bed  elevations  are  presented  at  x  -  0  (gate  location),  x  -  0.25  m  (up¬ 
stream),  and  x  -  0.25  m  (downstream).  It  is  seen  that  the  surface  elevation  decreases  upstream,  increases  downstream,  and 
remains  about  the  same  at  the  gate.  All  three  locations  experience  flushing  and  the  bed  elevation  decreases.  The  figure  shows 
that  there  is  reasonable  agreement  between  the  simulation  and  the  experimental  data  for  both  surface  and  bed  elevations.  A 
large  difference  between  simulation  and  experiment  occurs  mainly  at  the  initial  time  at  the  gate,  which  is  probably  attrib¬ 
utable  to  three-dimensional  effects  that  the  model  does  not  capture.  The  degree  of  agreement  between  numerical  simulation 
and  experimental  measurement  is  about  same  as  that  in  work  of  Abderrezzak  et  al.  [35]. 
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Fig.  7.  Solution  of  morphology  for  the  wave-driven  dune  problem 
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The  second  application  is  a  short  wave-driven  flow  over  a  sand  dune.  The  initial  conditions  for  the  wave  action,  the  veloc¬ 


ity  field,  and  the  bottom  shape  are,  respectively 

N  0.  (23a) 

U  0 ,  V  0  H  2  exp{  0.01  (x2  fy2)}.  (23b) 

Hb  exp{-0  01  (x2  t  y2)}.  (23c) 

The  upstream  boundary  condition  is 

H  2  m  x  =15  m.  (23d) 

The  wave  action  is  introduced  as  a  source  term  in  Eq.  ( 1 )  at  the  upstream  end 

S  =  0.15(1  -  tanh(20  4  x-O.Olt)),  (23e) 


which  generates  waves  starting  at  the  upstream  end  and  propagating  in  the  x-direction:  n  «  0.025,  af  «  1 ,  and  a  short  wave  is 
considered.  In  the  simulation,  Ax,  Ay  -  0.4  m,  Arr  -  0. 1  1  / s,  A 0  **  0.76  radians,  CFL  -  0.5,  and  von  -  0.2  Extrapolation  is  used 
to  update  solutions  at  the  boundaries. 

The  instantaneous  solutions  at  three  different  moments  are  given  in  Figs.  5-7.  As  shown  in  Fig.  7,  initially  (r «  427  s),  the 
shape  of  the  sand  dune  is  changing  into  a  triangle,  and  low-elevation  regions  are  forming  at  the  upstream  region,  lateral  sides 
and  immediately  downstream  of  the  sand  dune.  At  the  same  time,  the  sand  dune  is  propagating  downstream,  and  its  height 
gradually  decreases.  In  this  process,  the  sand  dune  becomes  asteroidal  (r *  500  s),  which  is  consistent  with  other  simulations 
(e.g.,  [16]).  The  low-elevation  regions  enlarge  and  merge,  and  the  elevations  at  the  upstream  region  and  lateral  sides  of  the 
sand  dune  become  the  lowest  (r  -  500  s,  568  s).  Finally,  the  sand  dune  is  removed  by  the  flow  and  the  final  seabed  elevation 
is  lower  than  the  initial  value  throughout  the  model  domain.  Fig.  6  shows  that  surface  elevation  reaches  its  highest  value  in 
front  of  the  sand  dune  but  forms  a  dip  immediately  behind  it,  and  it  gradually  decreases  in  all  regions.  Fig.  5  indicates  that 
wave  action  takes  larger  values  at  the  upstream  end,  and  it  propagates  downstream.  It  can  also  be  seen  from  Figs.  5  and  6 
that,  during  the  evolution  of  the  sand  dune,  both  the  wave  field  and  water  surface  elevation  evolve  congruently  with  the 
sand  dune.  Fig.  5  reveals  clear  traces  of  bed  morphology,  which  indicate  its  strong  effect  on  the  wave  field. 

5.  Concluding  remarks 

This  study  presents  a  framework  for  synchromcally  coupling  wave,  current,  sediment  transport,  and  bed  morphology  as 
well  as  preliminary  simulation  results.  The  numerical  algorithms  adopted  in  the  paper  lack  the  high-resolution  capability 
required  for  rapidly  evolving  flow  like  surges,  and  they  could  lead  to  spurious  oscillations.  Besides,  the  algorithms  have  se¬ 
vere  restrictions  on  time  step  because  of  their  explicit  forms.  Alternative  high-resolution  schemes  should  be  developed  for 
the  coupled  system.  More  comprehensive  validation  tests  with  experimental  and  field  data  are  needed  to  test  the  capability 
of  the  proposed  coupled  framework.  Furthermore,  in  order  to  provide  guidance  for  coupling  actual  models  (e.g.,  SWAN)  in 
practical  simulations,  it  would  be  useful  to  study  strategies  for  coupling  waves,  currents,  sediment  transport,  and  morphol¬ 
ogy,  that  will  permit  different  models  to  use  with  different  time  steps  and  grid  spacing. 
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